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• ' Quest for a new form of matter inside compact stars compels us to examine the thermody- 

f~^ ■ namical properties of the phase transitions. We closely consider the first-order phase transitions 

and the phase equilibrium on the basis of the Gibbs conditions, taking the liquid-gas phase 
transition in asymmetric nuclear matter as an example. Characteristic features of the mixed 
phase are figured out by solving the coupled equations for mean-fields and densities of con- 
stituent particles self-consistently within the Thomas- Fermi approximation. The mixed phase 
is inhomogeneous matter composed of two phases in equilibrium; it takes a crystalline structure 
with a unit of various geometrical shapes, inside of which one phase with a characteristic shape, 
called "pasta" , is embedded in another phase by some volume fraction. This framework enables 
us to properly take into account the Coulomb interaction and the interface energy, and thereby 
sometimes we see the mechanical instability of the geometric structures of the mixed phase. 
Thermal effect on the liquid-gas phase transition is also elucidated. 

Similarly hadron-quark deconfinement transition is studied in hyperonic matter, where the 
neutrino-trapping effect as well as the thermal effect is discussed in relation to the properties of 
the mixed phase. Specific features of the mixed phase are elucidated and the equation of state 
is presented. 
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9.1 Introduction 

Equation of state (EOS) for matter inside compact stars is sometimes characterized by the 
phase transitions. There have been expected various first-order phase transitions in nuclear 
matter, such as hquid-gas transition in the low densities, meson condensations and hadron-quark 
deconfinement transitions at high-density region [U EJ [31 [3j . If such phase transitions occurs 
inside compact stars, they affect EOS as well as thermal, dynamical or magnetic properties 
O El [71 E] • There has been a theoretical development and recent studies have shown that the 
inhomogeneous states called "pasta" may emerge as a mixed phase during the phase transitions, 
based on the Gibbs conditions [9]. Historically, the mixed phases have been mostly treated very 
naively by applying the Maxwell construction to get the EOS in thermodynamic equilibrium. 
Ravenhall et al. discussed the possible emergence of the pasta structures at low-density nuclear 
matter, since uniform nuclear matter is not stable at low densities [10]. Their pasta phase 
may be regarded as a mixed phase following liquid-gas phase transition. In 1990 Glendenning 
has pointed out that the phase equilibrium in multi-component system or in the case of plural 
chemical potential corresponding to the conserved quantities such as in neutron star matter 
must be carefully treated based on the more general Gibbs conditions [3l [11] . He emphasized 
that the phase transitions should be treated in a proper way in neutron-star matter, which is 
specified by two conserved quantities, electromagnetic charge as well as baryon number. In 
particular the local charge-neutrality condition used in the case of the Maxwell construction 
should be replaced by the global one in the mixed phase. Consequently the pressure becomes 
no more constant as in the Maxwell construction, but still increases in the mixed phase. He also 
demonstrated that the mixed phase is extended in the rather wide density region within a bulk 
calculation, where two sets of semi-infinite matter are considered in phase equilibrium, separated 
by a sharp boundary, and the pressure balance, chemical equilibrium and the global charge- 
neutrality conditions are imposed in the absence of the Coulomb interaction. It is important to 
see that particle fractions in each phase are no more the same to each other for given particle 
fractions as a whole. Nowadays we know that such phase transitions are rather common in other 
fields as well [HI [131 [HI [151 [161 [T7] a-nd we call them non-congruent phase transitions, following 
losilevskiy [17j . Generally any phase transition in a system with two or more conserved charges 
must be non-congruent. 

In this calculation the Coulomb interaction is discarded and we can by no means estimate 
how important it is, or the surface energy can not be taken into account. If we consider a more 
realistic situation, we must consider inhomogeneous matter with various geometric structures 
called pasta, and take into account the finite-size effects such as the Coulomb interaction and the 
surface tension. Heiselberg et al. have demonstrated that such finite-size effects are important 
for the mixed phase in the context of hadron-quark transition [18] . Treating the pasta structure 
as in ref. [10], they simply considered uniform particle densities in both phases. They showed 
a possibility of that the finite-size effects may largely restrict the pasta region. Subsequently, 
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Voskresensky et al. have shown that the rearrangement of particle densities and the screening 
effect for the Coulomb interaction becomes important, and the large surface tension gives the 
mechanical instability of the geometrical structures in the mixed phase [19]. Following this 
idea, we have studied the pasta structures brought by the kaon condensation and hadron-quark 
deconfinement transition [201 [2H [22] . More recently we have extended our framework to deal 
with the finite temperature case or the neutrino-trapping matter, which may be relevant to 
supernova explosions or neutron star mergers [23] . 

In this chapter we first present thermodynamic concepts of the non-congruent transition 
by considering the liquid-gas phase transition in asymmetric nuclear matter (section 19. 2p . In 
section 19.31 some results will be presented by a numerical calculation about the pasta structures 
in asymmetric nuclear matter. In section 19.41 we will present the results about hadron-quark 
deconfinement transition, where the effects of finite temperature or neutrino trapping are also 
discussed. Section [9.51 is devoted to summary and concluding remarks. 

9.2 Phase equilibria in thermodynamics 

9.2.1 Maxwell construction in the pure system 

Before beginning our discussion about the non-congruent transitions, let us briefly summarize 
more familiar case in the pure system, where the Maxwell rule can be applied to construct the 
equilibrium curve. Since there is only a single component, it is trivially congruent. Typical 
isotherms for liquid-gas phase transition are sketched in Fig. 19.11 where we can see that the 
gas phase (low-density phase) and the liquid phase (high-density phase), both of which are 
thermodynamically stable due to dP/dV < are clearly separated by the thermodynamically 
unstable region dP/dV > at temperature T below the critical temperature Tc. Above Tc 
two phases are indistinguishable]. The Gibbs conditions for phase equilibrium read in this case 
as P'^ = P^,iJp = pL^ besides Tq = Tl in terms of pressure Pi, chemical potential /ij and 
temperature Tj in each phase i = G or L. We can then construct the EOS by the graphical 
method, the Maxwell construction (see Fig. 19. ip . 

For a given T and fixed number, points A and B denote the initial and final points of the 
phase coexistence, subjected to the equal-area condition, 



V{P)dP= / V{P)dP, (9.1) 

AEM JbDM 

which is equivalent with ^ = /x by way of the thermodynamic relation, dG = Ndji = VdP 
with the Gibbs free energy G. Alternatively the double tangent method can be applied for 
the Helmholtz free energy F = G — PV . Anyway we can draw the binodal curve (coexistence 
curve) by connecting the equilibrium points, A and B, by changing T. At the same time we can 



* Recently the supercritical fluid above the critical point has attracted much interest, especially in the context 
of various applications [24] . 
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Figure 9.1: Isotherms of the van der Waals fluid on 
the P — V plane. 



Figure 9.2: Saturation curve on the P- 
T plane. 



also draw the spinodal curve by connecting the limit points of the stable region D and E where 
dP/dV = by changing T. At the point E, where d'^P/dV'^ > 0, the phase is intrinsically stable 
against the shrinkage of the volume, while it is unstable for the expansion of the volume. The 
similar argument can apply for the point D by exchanging the volume variations. Consequently, 
the segments AE and BD show the metastable regions. 

We can draw the saturation curve, which is defined by the equilibrium pressure as a function 
of T, on the pressure-temperature plane, as in Fig. 19.21 Then two phases are separated by a 
single curve terminated at the critical temperature T^, which means P is constant in the mixed 
phase. In the following section we shall see that it spreads to involve a finite mixed-phase region 
in the general case with multi-component, since pressure still increases in the mixed phase. More 
interestingly the mixed phase emerges with peculiar geometrical structures. 



9.2.2 Non-congruent transition and structured mixed phase within the bulk 
calculation 

First we briefly summarize the basic concepts of the non-congruent transition in a clear fashion by 
the bulk calculation without the Coulomb interaction or surface energy. We consider the liquid- 
gas transition at subnuclear densities in asymmetric nuclear matter ( ASNM) as an example [5S] . 
The proton-number fraction Yp = Np/N then specifies ASNM. 

It is well-known that EOS of nuclear matter exhibits a thermodynamically forbidden region, 
which separates "liquid" and "gas" phases at subnuclear densities. This phase transition looks 
like the first-order, so that we must treat the phase-equilibrium in the mixed phase. Since 
there are two conserved charges, electromagnetic charge and baryon number in ASNM, we must 
discuss the phase transition in the binary system. Different from the pure component system, 
thermodynamic behavior of the mixed phase becomes complicated in this case [24^ 126] : two 
coexisting phases can vary the chemical composition with keeping the total charges. 
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Thermodynamical potential for ASNM is defined as 
0(r, V, {fii}i=n,p) = -/3"Hn TV exp 
Then particle number and pressure are given as 



-/3 i? - J] fJ^iNi 



i=n,p 



(9.2) 



NdT,V,{H}n,p) = -{^) 



"^"'"'^'^^^-^^-(^)..M.; 



T,V,fij^i 



(9.3) 



The Helmholtz free energy can be expressed as 



FiV,T,Np,Nn) = VTiT,{piU,p) 



(9.4) 



i=n,p 



with the particle-number density pi , 



Ni = Vpi{T,{pj}n,p) 



(9.5) 



for uniform matter. The Gibbs-Duhem relation gives il = —PV. 

Usually EOS (P vs V) renders the Van der Waals type and thermodynamically unstable 
region (dP/dV < 0) appears at subnuclear density, which signals phase transition (called the 
"liquid-gas" phase transition in nuclear matter) [25j. It looks like first-order phase transition, so 
that there appears the mixed phase composing of two phases in equilibrium. Since there are two 
conserved charges, baryon number {N^ = Np + Nn = V ps) and isospin number {I^ = Np — N^ = 
Vpi), we must carefully treat the phase equilibrium in the mixed phase of this binary system^. 

In the following we assume the equilibrium between two phases (liquid(L) and gas(G)) in 
the mixed phase. Then the Gibbs conditions read 
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n,p). 



(9.6) 



The phase diagram is then given in the T-P-Yp space: phase-separation boundary is then 
given as a two dimensional surface. Note that the familiar Maxwell construction cannot be 
applied in this binary system, since it works only for single component system. When we draw 
the pressure surface as a function of /ip and /i„ for given Yp and T. we can graphically see a 
folded surface and the intersection gives the coexisting curve. The mixed phase develops along 
the coexisting curve and obviously satisfies the Gibbs conditions. As a result, pressure is no more 
constant in the mixed phase. So this is not first-order phase transition, but second-order one in 



In neutron-star matter electric charge is conserved instead of isospin at low temperature. Moreover, lepton 
number is also conserved in some situations as in supernovae, which makes a ternary system [36) . 



9.2.2 Non-congruent transition and structured mixed phase within the bulk calculation 



the sense of Ehrenfest |29j . However, since it never accompanies the change of any symmetry, 
it is first-order phase transition in the sense of Landau [26j . We shall see that the mixed phase 
consists of two distinct phases in equilibrium and there exist a metastable state and latent heat, 
which should characterize the first-order phase transition. 

When we plot the beginning and terminating pressures in the P-T plane by changing T for 
a given Y^, we have an extended coexisting domain. The phase-separation boundary is called 
hinodal curve. Following losilevskiy, we call it non-congruent phase transition. Note that pi 
as well as pB takes the different value in each phase, keeping the volume average fixed. It is 
instructive to seek the coexisting region by daring to apply the Maxwell construction, as in the 
van der Waals fluid (Fig. 19. ip . One may see only a curve on the same plane, which is called 
congruent phase transition and the water- vapor phase transition is a typical example [TTj . When 
we apply the Maxwell construction to the ASNM, the first two conditions in (j9.4p are obviously 
satisfied. However, the chemical equilibrium is partially violated since the proton fraction is 
forced to be the same in both phases: actually the baryon chemical potential, pB = Pp + Pn, 
results in being equal at the end points, while the isospin chemical potential, pj = pp — pn, is 
different. 

Consider the manifold J-" in the parameter space spanned by {pi}n,p with T fixed. If there 
is a common tangent plane connecting two distinct configurations with {p^}n,p in G and {p'-}n,p 
in L, Gibbs' criteria are satisfied. Then the mixed phase can be constructed by combining two 
phases with volume fraction A = V /V{= V + V ),0 < A < 1. Total charge densities are 
given by pi = Xp'l -|- (1 — A)/9^ and Helmholtz free energy can be written as 



HT, {p^]n,p) = A7-(r, {pnn,p) + (1 " X)HT, {Pi}n, 



(9.7) 
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Figure 9.3: Binodal curve for ASNM at T = lOMeV (see section [9.3. 2p . The filled circle denotes 
the critical point. 
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The binodal points for given Yp are then obtained by considering the two hmits, A = 0, 1 |j . 
From Eq. ([931) , 

P^(;U,)=P^(/i,), (9-8) 

and 

for A = 0, and 

QpL 

-^ = Pi (9.10) 

for A = 1. Three equations in (|9.4|) or Eqs. (j9.8|) . (19.91) . (19.101) determine the unknown values of 
chemical potentials /Xj and the baryon densities for the binodal points. Finally the binodal curve 
is obtained by connecting the binodal points by changing the mixture Y^ (see Fig. 19. 3p . Note 
that P^{fj.i) in the former case and P'^{^i) in the latter case never correspond to the A = and 
A = 1, respectively. 

9.2.3 Stability criteria 

Spinodal instability has been often discussed in ASNM in relation to the nuclear multi-fragmentation 
in heavy-ion collisions or in the crust region of neutron stars [SSI [301 ED 132 E31 E5] • Since it 
should be also related to the formation of the pasta structures, we briefly summarize it. The sta- 
bility criterion implies the convexity of the free energy density J-" of pj : Hessian or the curvature 

matrix 

/ d'^F/dpl d^T/dpndpp\ _ f dfin/dpn dpn/dpp \ , . 



d^Tjdppdpn d^Tjdpl ) \ dpp/dpn dpp/dpp 
has only positive eigenvalueqj, 

tr(C) > and det(C) > 0. (9.12) 

Since both eigenvalues never become negative in ASNM at subnuclear densities, one negative 
eigenvalue or equivalently det(C) < defines the spinodal region. The condition det(C) < is 
further equivalent with 

^^^ <0 or ft^l <0, (9.13) 



.dpB/TXp \9Yp/Tp 

where the first inequality is called "mechanical" instability and the second one "diffusive (chemical)" 
instability [31] . The mechanical instability is the familiar one and comes from the negative com- 
pressibility, and the spinodal condition is given only by this inequality in symmetric nuclear 
matter (SNM). However, the chemical instability first happens in ASNM, whereas curvature 
always remains positive and no spinodal instability is predicted in pure neutron matter [32]. So 



Note that there is an exception; we shall see a possibility of the existence of Amax in the special case, called 
retrograde condensation. 

'' In neutron-star matter, which is charge-neutral in the presence of electrons, the curvature matrix is modified 
by adding djj.e/dpe to dfj,p/dpp [501155] . 
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there should be a lower bound for Yp, below which ASNM is always stable. The Gibbs free energy 
G is depicted in Fig. 19.41 whose derivative with respect to proton number gives proton chemical 
potential, {dG/dNp)T^p = fJ-p- In Fig. 19.51 chemical potential isobar is shown for protons. In 
the region between the binodal points -61^2 the uniform matter is unstable and separated to two 
phases with different Yp, (0 < 1^ < 0.5) |f| ; the spinodal region is specified by the spinodal points 
5i^2 and two regions between Bi and 5i and B2 and 5*2 are the metastable regions. As ASNM 
enters the spinodal region, a unique mechanism of phase separation works and small fluctuations 
of Yp grow to produce a non-uniform random pattern. This phenomenon is known as spinodal 
decomposition [21]. In the metastable regions ASNM is stable for small fluctuations of Yp, but 
more likely becomes inhomogeneous by way of nucleation and growth. These features should be 
peculiar to the first-order phase transitions. 




^ip 




Figure 9.4: Schematic view of the 
Gibbs free energy G at fixed tempera- 
ture as a function of Yp. B12 denote 
the binodal points and 51^2 the spin- 
odal points where [d'^G/dYp) = 
NB{dfip/dYp)T,P = 0. 



Figure 9.5: Proton chemical-potential isobar at fixed 
temperature as a function of Yp. The segment 
between Si^2 denotes the spinodal region, where 
{d^p/dYp)rp p < 0, and the segments between Bi^2 
and 5i_2 correspond to the metastable regions. 



9.2.4 Finite-size effects and "Pasta" structures 

We have only discussed the binary system to demonstrate the characteristic features of the 
non-congruent transition. Since there are usually two conserved quantities, baryon number and 
electromagnetic charge in the cold catalyzed matter, many cases are treated as the binary systems 
inside compact stars. Extension of the discussion to multi-component system is straightforward 
but very complicated. In the neutrino trapping case as in supernova explosions, we need another 
conserved quantity, lepton number. Thus we must consider the ternary matter. In the recent 
paper this problem has been generally addressed and classification of various situations has been 
presented [36] . 



*The region 0.5 < i^ < 1 is redundant due to symmetry of interclianging protons and neutrons. 
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Figure 9.6: Schematic picture of "pasta" structures. When red phase appears in the blue phase, 
the red phase form droplets in the sea of blue. With increase of the volume fraction of red phase 
in the blue phase, the shape of the red phase changes from droplet to rod, slab, tube, and to 
bubble before pure red phase appears. 

We have presented general features of the non-congruent transition in a clear fashion with 
a help of the bulk calculation, discarding the Coulomb interaction and geometrical structures 
which are called the finite-size effects. However, the bulk calculation is insufficient to describe 
the essential aspects of the non-congruent transition. We can see the emergence of various geo- 
metrical shapes, called pasta (see Fig. I9.6p . by the cooperative effect of the Coulomb interaction 
and the surface energies, and the importance of the charge screening effect which is a typical 
many-body effect and makes pasta structures unstable, called mechanical instability. 
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Figure 9.7: Wigner-Seitz cell for the droplet case (left panel) and crystalline structure composed 
of droplets (right panel). 



A rough estimate about the finite-size effects can be done as a straight extension of the bulk 
calculation. Since the crystalline structure appears in this case we use the Wigner-Seitz cell 
approximation (see Fig. 19. 7p . Consider a simple model with a droplet composed of the phase I 
embedded in the phase II in a Wigner-Seitz cell with the volume Vw = AttR'^/S at T = 0: particle 
density is assumed to be uniform in each phase, charge neutrality is achieved by electrons and 
their interface is represented by the sharp surface with a surface tension parameter ts- For a 
given volume fraction A = (R/Rw)^, the total energy then consists of three parts. 



E = Ey + Ec + Es, 



(9.14) 
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where Ey is the volume energy computed with the strong-interaction Hamiltonian, and Eq and 
Es are the Coulomb energy and the surface energy, respectively. The energy density Ey/Vw 
does not depend on R, while the Coulomb-energy density renders 



Ec/Vw = ^~[^(pih - pihfR^, 



and the surface-energy density 



Es/Vw = A 



375 
R ' 



(9.15) 



(9.16) 



In Fig. 19. 81 we sketch the R dependence of these energy densities. The optimal value of R is then 

given by the minimum condition, 

d{E/Vw) 



dR 



0. 



(9.17) 



Using Eqs. (|9.14p . (|9.15|) . (|9.16|) . we can see the well-known relation Es = 2Ec holds at the 



minimum. 




Figure 9.8: Energy densities as functions of the radius of droplets. Solid curves show the case 
without the screening effect, while dotted curves with screening. 



We shall see this simple consideration is too crude to describe the mixed phase in a realistic 
way. Once the Coulomb interaction is included, we must take into account the rearrangement 
of the particle densities at the same time, and in turn the screening effect for the Coulomb 
interaction. Consequently the Coulomb interaction becomes effectively short-range by the Debye 
mass, mj^ = Yli'^'^'^Plh/'^Pi- ^^ shown in Fig. 19.81 the Coulomb-energy density is no more an 
increasing function with respect to R, but is exponentially reduced at large R. The minimum 
point is shifted to larger R as ts increases, which gives rise to a mechanical instability of the 
droplet above the critical value of ts- The mechanical instability of the pasta structure also 
occurs in other geometrical shapes, which means that both R and Rw go to infinity with their 
ratio fixed. Thus we recover the picture of the phase equilibrium of two bulk matters, where 
the surface energy is irrelevant and the Coulomb energy vanishes due to the achievement of the 
local charge-neutrality. These features indicate that the EOS resembles the one given by the 
Maxwell construction. 
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The rearrangement of the particle densities also induces a change in the volume energy Ey , 
which is composed of kinetic energy and strong-interaction energy. In the above calculation Ey 
itself is intrinsically large but independent on R, since we have assumed uniform densities p- ' 
in each phase. Once the Coulomb interaction is properly taken into account, particle densities 
are no more uniform to produce a R dependence in Ey- The difference in Ey is called the 
correlation energy in ref. [19], defined by 

Ecorr = Ey{pi{v)) - Ey{pl'^^). (9.18) 

The detailed analysis showed that -Bcorr gives negative values and behaves R^^ for large R. 
Sometimes .Bcorr gives a significant contribution as later seen in hadron-quark deconfinement 
transition (section 19. 4p . 

9.3 Liquid-gas transition in asymmetric nuclear matter 

The properties of inhomogeneous low-density matter was already studied in the early 70th by 
Baym et al. [31] • They have presented a concept of structured mixed phase which includes 
nuclear droplet in the neutron gas as well as its reversal structure (droplet of gas in the nuclear 
sea). 

About a decade later, Ravenhall et al. [lOj have considered not only the three-dimensional 
structures but also those of one and two dimension. They have introduced a dimensionality d 
with a fractional number as well as other parameters such as volume fraction, density in the 
nuclear phase, and so on. With given average density and proton fraction, they looked for the 
parameters which give a minimum free energy. With increasing the average density, the structure 
of nucleus changes from droplet to rod, slab, tube, bubble, and finally uniform. Afterward, such 
a series of structures was called "pasta" (see Fig. 19. 6p . The corresponding d changes from 3 for 
droplet to 2 for rod, 1 for slab, 2 for tube, and 3 for bubble. But d changes smoothly, which 
suggests the existence of intermediate shapes besides the typical pasta structures. 

The next year, Hashimoto et al. [35] used a liquid-drop model with the Wigner-Seitz (WS) 
cell approximation, and showed that the stable nuclear shape changes as in the same way 
as has shown by the above Ref. pJDj. After that there appeared studies of pasta structures 
by several groups using different frameworks which are, however, basically mean-field models 

[MiiiQiiiiiJiaiiiiiii]. 

Besides the mean-field approaches, molecular dynamics (MD) simulations \Ab\ H6] were ap- 
plied to the inhomogeneous low-density matter [TJ HU HTJ HH] . The advantage of MD simulation 
is that one needs to assume no structure, or capability in simulating the dynamical processes [7]. 
Some mean-field approaches are also possible without the assumption of the structure [iHl [50] . 



9.3.1 Framework based on relativistic mean-field theory 



12 



9.3.1 Framework based on relativistic mean-field theory 

Here we use the relativistic mean-field (RMF) model to describe asymmetric nuclear matter. 
Exploiting the idea of the density functional theory (DFT) within RMF model [51], we can 
calculate the particle-density profiles as well as thermodynamic quantities. The Lagrangian 
density of the system is written as 






C-N + C,M + Ce 
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where -F^i^ stands for d^Fi, — dyF^. By the Euler-Lagrange equation 



d„ 



dC 



d{d. 



dC 



0, 



cr,W/„i?/„V;„^) 



(9.19) 
(9.20) 



(9.21) 

(9.22) 
(9.23) 
(9.24) 



(9.25) 



we get the coupled equations of motion for the scalar (cr), isoscalar- vector (w), isovector- vector 
(p), electromagnetic and nucleon fields. Since the spatial components of vector fields become 
zero due to the isotropy of the system, the meson mean- fields (cr, a;o,i?o) and electric field (Vq) 
obey the following equations of motion 



- V^o-(r) + m^o-(r) 

-V^tJo(r) +m^tJo(r) 

-V2flo(r) + m2i?o(r) 

VVo(r) 



dU 



+ g.^(pW(r)+pW(r)) 



d(T 

gcvN{pp{r) +Pn(r)) 

9pN{pp{r) - p„(r)) 

47re^Pch(r) (Pch(r) = Pp(r) - Pe(r)). 



(9.26) 

(9.27) 
(9.28) 
(9.29) 



For nucleons and electrons we employ the Thomas-Fermi approximation at finite temperature 
T and the distribution function / of proton, neutron and electron can be written as 



/i=n,p(r;p,A*j) = (1 + exp 

/e(r;p,/ie) = (l+exp[(p-(^e-V^o(r)))/T])-i. 



^Jp^+mUrf-lyi{r))/T 



-1 



(9.30) 

(9.31) 
(9.32) 
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For simplicity we neglect contributions of anti-particles, which is legitimate for high chemical 
potentials /ij and low T. Using above distribution functions fa, {a = n,p,e), relations between 
chemical potentials fia and densities pa are written as 

Pe(r) = 2/ -^;^/e(r;p,^e), (9.33) 

pi=p,n{r) = 2/ -— 3/i(r;p,/ii), (9.34) 

Pn = i^n{r) +gu>N^o{r) - gpNRoir), (9.35) 

Pp = i/p(r) +sr^Ara;o(r) +5rpAri?o(r) - Vb(r), (9.36) 

We demand that the chemical potentials pa (a = p, n, e) become independent of the position 
r. Note that the condition Pp = Pn — Pe among chemical potentials is applied only in the case 
of beta-equilibrium. Here we'd like to make a comment about the inclusion of the Coulomb 
potential. The electron chemical potential as well as the proton chemical potential is the gauge- 
variant quantity and only the charge density is gauge-invariant. Hence, the electron density is 
not uniform in the presence of the Coulomb potential even when the condition pe =const. holds 
everywhere |19j . 

Accordingly we list in addition some necessary quantities, the scalar density p^^' used in 
Eq. (j9.26p . energy density e and the entropy density s below: 

pUpJr) = !I!^Hs{r;u,,mUr)) (9.37) 

i=n,p 

+\ (Vwo(r))' + ujo{r){pp + p^) - ^uj^r) 

1 771 

+- (Vi?o(r))' + 5p7vi2o(r)(pp(r) - p„(r)) - -^i?^(r) 

^ ^9pN 

^\ (Va(r))^ + ^a\.) + ^<.3(r) + |,4(r) 



1 /"°° 
H — 2/ /^P/e(r;p,/Ue) +ec(r) 
"^ JO 

/■°° Sp 

Sa=n,pA'^\Pa) = 2/ — — 3 [/„(r; p, /Xa) log (/a(r; p, ^a)) 



+ (1 - /a(r; p, Pa)) log (1 - /a(r; p, pa))\ , (9.38) 

where £c is the contribution from the Coulomb interactioro, 

ecij) = ^ / ci3,. /^ch(r)peh(rO ^ ^^^^^ 

^ Jv 1"^ ~ r I 



® We here discarded the exchange term, which is evaluated in the local-density approximation as e™ = 
-3e74(3/7r)i/Vch(r)^/^ M- 
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Table 9.1: Parameter set used in RMF in our calculation. With these parameters, the saturation 
property of nuclear matter is reproduced: minimum energy per baryon —16.3 MeV at p = 0.153 
fm~^ with the incompressibility K = 240 MeV. 
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9pN 
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6.3935 


8.7207 


4.2696 


0.008659 


-0.002421 


400 


783 


769 



and the functions Hn-, (n = 3, 5) are defined by 



Hn{r;ui,m) = / — — — y^ (/^(r; p,^) - /i(r; p, -/x)) 



^" "^^ ;fi{r■,P,^^)■ (9.40) 



(p2 + m2)i/2- 

To solve the coupled equations of motion (J9.26p - (j9.36|) . we employ the Wigner-Seitz approx- 
imation, where whole the space is divided into equivalent cells with a geometrical symmetry in 
one, two or three dimension. The cell is divided into grid points and the coupled equations of 
motion are solved numerically. The size of the WS cell is optimized so that the Helmholtz free 
energy. 



Fws 



f d^rle-T Y^ sa] (9.41) 

becomes minimum, where V\vs denotes the volume of the WS cell. In the following we consider 
two cases: one is a given proton fraction case, Yp = jy Ppi^)/ jy iPpi^) + Pn{^)) and the 
other is /3-equilibrium case with additional condition, pn = Pp + Pe The value of pe is adjusted 
at any time step to maintain the global charge neutrality, Jy d^rpchi^) = 0: we decrease pe 
when the total charge in a cell is positive and vice versa. All the above relaxation procedures 
are carried out simultaneously. The parameters used in our calculation are fitted to give the 
saturation properties of nuclear matter at T = and properties of finite nuclei in the ground 
states |53j, which are listed in Table I. 

9.3.2 Bulk Gibbs calculation 

In this section we discuss about bulk properties of nuclear matter without taking into account 
of electrons and the finite-size effects. We give some isothermal curves in the Ps-Yp plane and 
isolepthal curves in the Pb-T plane by the bulk Gibbs calculation, where Pb denotes the baryon 
partial pressure. Dropping the electric field V and assuming uniform densities Pi,{i = n,p) 
as well as the meson mean-fields, the thermodynamic potential density and the energy density 
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(9.43) 



where we subtracted a rather trivial electron contribution. The entropy density is given by 
Eq. (j9.42p or follows from the Gibbs-Duhem relation, 



SB = {Pb + ^b - X^ PiPi)/T, 



i=n,p 



SO that we have the Helmholtz free energy density, 

^B = ^B — Tsb- 



(9.44) 



(9.45) 




Figure 9.9: Temperature-dependence of phase coexistence curves in Yp-Ps plane for SNM 
Yp = 0.5. Filled circles show the critical points. 
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Figure 9.10: Same as Fig. [9J] for T = 10 MeV. 



In Fig. 19.91 binodal curves with several temperatures are presented in the Yp-Ps plane. For 
a given Yp and temperature T, e.g. Yp = 0.3 and T = 10 MeV (Fig. [gTTOl ASNM is in the 
gas phase when pressure is small enough. As pressure increases the system crosses the binodal 
curve, where the phase separation occurs and mixed phase begins (A = 0). The gas and liquid 
phases coexist in the mixed phase and the proton fraction as well as baryon number density 
in each phase are different. As pressure further increases, the volume fraction of each phase 
A get larger, and eventually total system becomes the pure liquid phase (A = 1). The critical 
point is given by the condition, dPs/dYplcp = 0, since the two phases cannot be distinguished 
there, Pb{T,YS^^) = Pb{T,YS'^ + 5Yp). On the other hand, if we dare to apply the Maxwell 
construction (forced congruent transition), we have the congruent curve, where the critical point 
is the terminal point. 

Note that the phase transition is congruent for symmetric nuclear matter (1^ = 0.5). This 
is because the free energy should be smallest at 1^ = 0.5 for every density. In this case the 
trajectory of the system and the coexistence curve meet at a single point on the Yp-Ps plane 
and the mixed phase of this system consists of gas and liquid with the same value of 1^ = 0.5. 
Accordingly dPs/dYp = at this point (Gibbs-Konovalov theorem 



In the case of 1^ = 0.1, the system starts as a gas phase in the panel at the far right of 
Fig. 19.101 With increase of density, the pressure goes up till Pb ~ 0.08 MeV • fm~^. Then 
the LG mixed phase appears as in the central panel. With further increase of density, the 
pressure goes up until Pb ~ 0.18 MeV • fm~^. At this point Yp of the gas phase becomes that 
of the total system, 0.1, which means that the system becomes a gas again. However such 
"supercritical" gas has very different properties from the usual gas at low pressure. Such a 
transition, (gas) — > (mixed phase) -^ (supercritical gas), is peculiar to systems with two or more 
chemical components and called "retrograde condensation" |26) . 

This situation is more clearly seen in Fig. 19.111 , where the binodal curves are shown in the 
Pb-T plane for given Yp. 

The critical point (i^cp^^^cp) is located on the binodal curve, but it is not the special point 
on the curve. The supercritical region (T > Tcp) also appears in this plane for 1^ = 0.1 and 
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Figure 9.11: Phase coexistence curve on T-Pb plane. Gray area shows the region of LG mixed 
phase. Dots on the curves indicate the critical points. The panel at the far left exhibits the 
congruent curve (see Fig. I9.2p . while others the non-congruent ones. 

0.2. 

There is no area inside the coexistence curve for Yp = 0.5, since the system behaves congru- 
ently. For smaller 1^, on the other hand, the area become wider. The gas phase for Yp < 0.2 
includes dripped neutrons even at T = 0. therefore, the curve is not closed at T = 0. 



9.3.3 Pasta and EOS 

In this section we discuss the non-uniform nuclear matter by taking into account the finite-size 
effects. We include electrons to make the total system to be charge-neutral. Electrons are 
distributed quasi uniformly even in the mixed phase and the finite-size effects by electrons is 
small. We can see that pasta structures with various geometrical shapes appear in the mixed 
phase. The properties of asymmetric nuclear matter with pasta structures will be elucidated by 
comparing them with those given by the bulk Gibbs calculation in the previous section. 

We show in Fig. 19. 12] some typical density profiles in the cell. The left and the middle panels 
show the cases of proton fraction Y^ = 0.1 and Yp = 0.5. Apparently, dense nuclear phase (liquid) 
and dilute nuclear /electron phase (gas) are separated in space and they form pasta structures 
depending on density. One should notice that coexisting two phases are non-congruent and 
have different components, i.e. nuclear matter and electron gas. Therefore the EOS of the whole 
system cannot be obtained by the Maxwell construction. Since electron density is almost uniform 
and independent of baryon density distribution, we can separately discuss the properties of the 
baryon partial system. 

In the case of 1^ = 0.5, as discussed in Sec. 19.3.21 the local proton fraction stays almost 
constant. Therefore the system behaves like a system with single component. This means that 
one can apply the Maxwell construction to get the baryon partial pressure Pg as in the left 
panel of Fig. 19.131 uniform low-density matter with a negative partial pressure is not favored 
and the Maxwell construction gives Pb = for the mixed phase. By the finite-size effects, i.e. 
the Coulomb repulsion and the surface tension, the structured mixed phase becomes unstable 
in the density region just below po, and consequently uniform matter with a negative partial 
pressure is allowed. Note again that the total pressure including electrons is always positive 
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Figure 9.12: Examples of the density profiles in the Wigner-Seitz cells, with 1^ = 0.1 (left), 
Yp = 0.5T = (middle) and T = 10 MeV (right). 

even in this case, so that the system is thermodynamically stable. 

In the case of asymmetric nuclear matter, e.g. Yp = 0.1 in the left panel of Fig. 19.121 the 
proton fraction in the dilute and dense phases are different, especially for low Yp. Matter behaves 
as a system with multi chemical components and the Maxwell construction does not satisfy the 
Gibbs conditions. 



9.3.4 Thermal effects 

Let us discuss the thermal effects on the LG mixed phase of low density nuclear matter. By 
comparing the density profiles in the middle and right panels of Fig. 19.121 we easily notice that 
the dilute phase at finite temperature always contains baryons while they are absent in the 
dilute phase at zero temperature if 1^ ~ 0.5. This is due to the Fermi-Dirac distribution at 
finite temperature, where density as a function of chemical potential is always positive. 

We also notice that the size of the pasta structure is smaller in the case of finite temperature. 
This comes from a reduction of the surface tension between two phases at finite temperature 
since the difference of baryon densities between two phases get smaller. 

The isothermal EOS's (baryon partial pressure as a function of baryon number density) 
of symmetric nuclear matter at various temperatures are shown in the left panel of Fig. 19.131 
Dotted and thick solid curves show the cases of uniform matter, while thin solid curves are 
the cases where non- uniform pasta structures are present. As shown in the right panel of Fig. 
19.131 pasta structures appear at finite temperatures as well as the case of T = 0. But there 
appears uniform matter (gas phase) at the lowest-density region |27tl28j since the baryon partial 
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Figure 9.13: Baryon partial pressure as a function of density pB for symmetric nuclear matter 
Yp = 0.5, which is given by subtracting electrons pressure from the total one. 

pressure of uniform matter has a positive gradient against density. On the other hand, the 
uniform matter is unstable where the pressure gradient is negative even if the pressure itself is 
positive. At T = 14 MeV, we obtain no pasta structure since the baryon partial pressure of 
uniform matter becomes a monotonic function of density above this temperature. Study of the 
instability of uniform matter and the appearance of the pasta structures in connection with the 
spinodal region is in progress. 

In Sec. 19.3.21 we have shown a retrograde condensation at T = 10 MeV and 1^ = 0.1 by the 
bulk calculation. In the full calculation with pasta structures, however, there is no evidence to 
have such phenomenon so far. Probably it is washed away by the finite-size effects. 



9.3.5 Discontinuity in pressure 

If one applies bulk Gibbs conditions to the low-density nuclear matter, the pressure changes 
continuously with the density. Mechanically unstable region with negative gradient of pressure 
dP/dpB < never appears. However, with inclusion of the finite-size effects for inhomogeneous 
nuclear pasta, some part of density region becomes unstable for the mixed phase. In this density 
region, uniform matter with negative baryon pressure is allowed (see Fig. I9.13p . The baryon 
pressure of the mixed phase also can be negative. However, the total pressure including that of 
electron and its gradient are always positive. 

Another thing to note is that there are discontinuities in pressure at the densities where the 
geometrical structure changes. The free-energy density as a function of density is continuous 
since it is the criterion of choosing the optimal structure. On the other hand the pressure Pb = 
^j PiPi — e is not necessarily continuous since the chemical potential pi should be dependent on 
the structure. So far we employ the Wigner-Seitz cell approximation for matter with the pasta 
structures. This is rather a strong constraint on the structure. In the work of Ravenhall et al. 
|10j dimension parameter is optimized to get the ground state structure. It behaves smoothly at 
the transient region except for a jump between uniform and bubble. If we could consider more 
than the typical pasta structures |50j . the discontinuity of pressure may diminish. 
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9.4 Hadron-quark deconfinement transition 

In this section we discuss the emergence of the pasta structure associated with the hadron-quark 
deconfinement transition ^i. 

Nowadays the hadron-quark deconfinement transition and its impfications has been exten- 
sively studied by various viewpoints: a great deal of effort has been paid for quest of signals 
in relativistic heavy-ion collisions [SB], compact-star phenomena [S] or remnants of primordial 
QCD transition [57]. Theoretically, numerical simulations based on the lattice gauge theory 
have been done to estimate the critical temperature and figure out the properties of the phase 
transition. However, for the present, any numerical method is useless for the deconfinement in 
compact stars. Instead, the studies by using the effective models of QCD such as the MIT bag 
model or the NJL model have been done. 

Although the order of the phase transition is not clear, theoretical arguments have suggested 
it should be of the first order in the low temperature and high-density regime |56j . So, let us 
assume here it is of the first order. Then we can expect the emergence of the structured mixed 
phase as in liquid-gas transition. Actually Glendenning first discussed it to demonstrate how the 
Maxwell construction is insufficient in neutron-star matter. Before the appearance of his paper, 
when the first-order phase transition occurs and new phase appears, EOS has been constructed 
by simply applying the Maxwell construction. Using the bulk calculation without the Coulomb 
interaction or the surface energy, called the finite-size effects, he could show a wide density 
region of the mixed phase for hadron-quark deconfinement transition, the pressure of which is 
by no means constant. 

Once the finite-size effects are properly taken into account, inhomogeneous matter called 
pasta emerges as is already shown in the previous section. Heiselberg et al. immediately 
pointed out the importance of the finite-size effects and demonstrated how the surface en- 
ergy parametrized by the surface-tension parameter makes the region of the structured mixed 
phase narrower [18] : they showed the disappearance of the structured mixed phase for a critical 
value of the surface tension. They considered only droplets for hadron-quark deconfinement 
transition, while more general pasta structures or the crystalline structures have been subse- 
quently discussed [55]. Voskresensky et al. discussed the charge screening effect on the stability 
of the structured mixed phase during the deconfinement transition, and showed its relevance 
for a given surface tension |19j . They emphasized that charge screening and rearrangement 
of particle densities become important for the description of the mixed phase and EOS. Using 
the effective model for hadron matter and the MIT bag model for quark matter, Endo et al. 
explicitly demonstrated these features [21j . 

The surface tension should be calculated self-consistently with models of hadron and quark 
matter. Unfortunately, there is a large ambiguity about its value, whereas a large value of 
100MeV-fm~^ might be more realistic [54] ■ 



^° We, hereafter, dismiss the possibihty of color superconductivity 
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9.4.1 Framework using Brueckner-Hartree-Fock calculation and the MIT bag 
model 

There are many models to describe the EOS with hyperons. However, nowadays, the potentials 
of the baryon-baryon interactions have been directly calculated by the great successes in the 
lattice QCD calculations [58]. Hence, it is very hot topic to calculate the EOS from the baryon- 
baryon interaction directly. Here, we show the EOS by the theoretical framework of the nonrela- 
tivistic Brueckner-Hartree-Fock approach |59j based on the microscopic nucleon-nucleon {NN), 
nucleon-hyperon (NY). The Brueckner-Hartree-Fock calculation is a reliable and well-controlled 
theoretical approach for the study of dense baryonic matter. Detailed procedure can be found 
inRefs. pi EQI EH [62] . 

We adopt the Argonne Vig potential [63j for A^A^ potentials, and semi-phenomenological 
Urbana UIX nucleonic three-body forces [M] and the Nijmegen soft-core NSC89 NY potentials 
[65| . Unfortunately, there are not reliable potentials for the hyperon-hyperon interaction now. 
Therefore, we neglect them here. We strongly hope that the lattice QCD calculations or the 
experiments, such as J-PARC, will reveal the these interactions in the future. 

Moreover, we discuss the hadron-quark mixed phase with the finite-size effects at finite tem- 
perature with/without trapped neutrinos. To take into account the effects of finite temperature, 
we adopt Frozen Correlations Approximation [MliSZliSH]. In this approximation, the correla- 
tions at finite temperature are assumed to be the same with the ones at zero temperature. It is 
found to be good accuracy at finite temperature by past studies [Ml [671 (SB]- Accordingly, we 
focus on the hadron-quark mixed phase at finite temperature and adopt this approximation. 

In the following, we briefly describe our framework. First we get the chemical potential /Xj 
from the number density rii, 



CO 



Pi = j^j Mp)^7Tp'dp (9.46) 

fiip) = 77 V7¥x-7^ (9-47) 

exp{(ei - Hi)/T} + 1 

where Si and fi{p) are the single-particle energy and the Fermi-Dirac distribution function, 
respectively, whereas subscript i denotes the particle species, i = n,p, A,S^. We set each 
degeneracy factor g = 2, and adopt each mass as m„ = nip = 939 MeV, m\ = 1115.7 MeV, and 
m^- = 1197.4 MeV. We note that Si includes the interaction energy Ui as well as the kinetic 
energy [661169]. 



e^ = ^m^+p^ + Ui. (9.48) 

Finally, we get the free-energy density J-" as following 



-^^ = E I (^/" \/^1+P^Mp) ^^P'dp + \ U^p^ I - TSH, 



(9.49) 
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where sh is the entropy density calculated from 

'H = -T^J^J {fi(p)^^MP) + (1 - Mp))H^ - Mp))} 4V^p. (9.50) 

The total pressure for uniform hadron phase is given by 

PH = Y.^^^Pi-J'H. (9.51) 

i 

Next, we show the quark matter based on the thermodynamic bag model. In this model, we 
get the number density ng, the pressure Pq, and the energy density eg. 



PQ = Y.J^l f<iip)Wdp, (9.52) 

^Q = T^jSasJ ^'nril+p^fq{p)AT,p''dp + B, (9.53) 

'Q = -Y.J^J^ {fg{p)lnfg{p) + {l-fgip))Hl-fg{p))}Wdp (9.54) 

^Q = eg- TsQ (9.55) 

for uniform quark matter, where fq{p) is the Fermi-Dirac distribution function of the quark q 
(= u,d,s), ruq its current mass, and B the energy-density difference between the perturbative 
vacuum and the true vacuum, the bag constant B set to be 100 MeV-fm~^. 

However, there are many uncertainties for B. Here, we also adopt the another type of B by 
the density dependent bag model. In this model, we set the vacuum energy B{pb) as 

B{pb) = Soo + {Bo - 5oo) exp[-/3^] 

Po 

with 5oo = 50 MeV-fm"^, Bq = 400 MeV-fm"^, and /? = 0.17. Here, we set the saturation 
density as po = 0.17 fm^'^. This approach has been proposed in Ref. [68], on the basis of 
experimental results obtained at CERN SPS on the formation of a quark-gluon plasma. The 
just mentioned choice of the parameters allows symmetric nuclear matter to be in the pure 
hadronic phase at low densities, and in the quark phase at large densities, while the transition 
density is taken as a parameter. 

For the mixed phase, we must numerically search the minimum of the free energy 

F= j dr^TH[n^] + I dr'^Tqiuq] + F, + F,^ + Ec + Eg (9.57) 

Jvh Jvq 

changing the volume ratio of quark matter to hadron matter, the geometrical types of structures, 
and the size of Wigner-Seitz cell Rw- Here, the free energies of electrons, Fg, and neutrinos, 
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Figure 9.14: Density profiles and Coulomb potential Vc within a 3D (quark droplet) Wigner- 
Seitz cell of the mixed phase at pB = 0.4 fm~^. The cell radius and the droplet radius are 
Rw = 26.7 fm and R = 17.3 fm, respectively. 



Fi,^ , are calculated by the use of the Fermi-Dirac distribution functions. We do not consider all 
kinds of anti-particles for simplicity. 

In the above equation, Es stands for the surface energy which comes from a sharp boundary 
between hadron and quark phases with a fixed surface tension rg. We demand the global charge 
neutrality and the chemical equilibrium as, 



IJ'U + fJ-e — Pue = /^d = /^s, 

Hp + He- fJ-Ue = fJ-n = PA = fJ-u + 2/Urf, 



(9.58) 



Although the surface tension of the hadron-quark interface is poorly known, some theoretical 
estimates based on the MIT bag model for strangelets [70] and lattice gauge simulations at finite 
temperature [71] suggest a range of ts ~ 10-100 MeV-fm~^. We will discuss the effects of surface 
tension on the pasta structures in subsection 1.4.4. 

Fig. 19. 141 shows an example of the resulting density profile in the droplet phase for pB = 0.4 
fm~^ at T = 0. One can see the non-uniform density distributions of particle species together 
with the non- vanishing Coulomb potential. The quark phase is negatively charged, so that d and 
s quarks are repelled to the phase boundary, while u quarks gather at the center. The protons 
in the hadron phase are attracted by the negatively charged quark phase, while the electrons 
are repelled. 

9.4.2 Hadron-quark mixed phase and compact stars 

Once obtained the EOS of matter, we can know the density profile (density as a function of 
radial distance from the center) and consequently the mass and the radius of compact stars. 
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Figure 9.15: Left panel shows the EOS of the mixed phase (thick dots) in comparison with pure 
hadron and quark phases (thin curves) at zero temperature. We also show, for comparison, the 
mixed phase by the Maxwell construction by thin solid line. Note that the Maxwell construction 
is incorrect in the present case, but it is still useful as an eye guideline. Right panel show mass- 
radius relation of calculated neutron stars. 



This is done by numerically solving the Tolman-Oppenheimer-Volkoff (TOV) equations pLj as 

d-p _ Gme (1 +p/e) (l + ATxr^p/m) 



dr 
dm 
dr 



1 - 2Gm/r 



—— = Anr e , 



(9.59) 
(9.60) 



being G the gravitational constant. Starting with a central mass density e(r = 0) = €c, one 
integrates out Eqs. (|9.59p and (|9.60p until the surface pressure becomes to zero. This gives the 
stellar radius R and its gravitational mass M = m[R). 

We show the relation between EOSs and structures of compact stars in Fig. 19.151 The left 
panel is the EOSs and the right panel shows the mass-radius relation of calculated neutron stars. 
As for the EOSs, we shows three types of quark-hadron phase transitions; e.g. the EOSs by the 
bulk Gibbs condition, the Maxwell construction, and full calculation considering the finite size 
effects. Here, the surface tension is set as ts = 40 MeV-fm~'^. We also show the uniform hadron 
and quark matters for comparison. Using these EOSs and the TOV equations, we can calculate 
the mass-radius relations shown in the right panel. It is well-known that the observed masses 
of neutron stars are distributed in a narrow band close to 1.4 times the solar mass Mq |72j . 
The maximum mass of neutron star is enough larger than the observed masses. With inclusion 
of hyperon degree of freedom, however, the softening of matter reduces the possible mass of 
neutron stars as low as 1.3 Mq, which contradicts the observation. The blue curves show the 
cases with hadron-quark mixed phase. The upper most one and the bottom are the results with 
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Figure 9.16: (a) Particle fractions of neutral matter with electrons, (b) The same quantity 
for charged matter without electrons, the low-density part of which corresponds to symmetric 
nuclear matter. 



the EOSs obtained by the Maxwell construction and the bulk Gibbs calculation, respectively. 
The result with inclusion of the full pasta structures with Tg = 40 MeV • fm~^ lies just below 
the Maxwell construction case. 

In any case with the hadron-quark mixed phase, the calculated maximum mass is slightly 
above 1.4Mq, which might be tolerable in confrontation with the observation. What is important 
is that the mixture of hyperons in nuclear matter softens the EOS, which causes the contradiction 
with the observation, i.e. the "maximum-mass problem". By considering the deconfinement 
transition to quark matter, the EOS gets stiffer at higher density region. This stiffness at high 
density sustains the neutron stars with M ~ IAMq. 

9.4.3 Hyperon suppression in the mixed phase 

Though we allow the freedom of hyperon mixture, hyperons don't actually come into the mixed 
phase. The suppression of hyperon mixture in the mixed phase is due to the fact that the 
hadron phase is positively charged. As shown in Fig. 19.161 hyperons (S~) appear in charge- 
neutral hadronic matter at low density to reduce the Fermi energies of electron and neutron. 
Without the charge-neutrality condition, on the other hand, there is symmetric nuclear matter 
at lower density and hyperons will be mixed there at higher density due to the large hyperon 
masses. Generally speaking hyperons are hardly mixed in the positively charged matter. Thus, 
the mixture of hyperons is suppressed in the mixed phase where the hadron phase is positively 
charged. 

9.4.4 Effects of surface tension 

Let us consider the role of the surface tension on the mixed phase. As already mentioned, 
the surface tension between hadron and quark phases is unknown. If one uses a smaller surface 
tension parameter rg, the energy gets lower and the density range of the mixed phase gets wider. 
The limit of r5 = leads to a bulk application of the Gibbs conditions without the Coulomb 
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Figure 9.17: Droplet radius (R) dependence of the energy per baryon for fixed baryon density 
Pb = 0.35 fm~^ and different surface tensions. The temperature T is zero for all cases. The 
quark volume fraction {R/Ry/Y is fixed for each curve. Dots on the curves show the local energy 
minima. The black line shows the energy of the Maxwell construction case. 



and surface effects. On the other hand, using a larger value of Tg, the geometrical structures 
increase in size and the EOS gets closer to that of the Maxwell construction case. Above a 
limiting value of ts ~ 65 MeV • fm~^ the structure of the mixed phase becomes mechanically 
unstable[19]: for a fixed volume fraction {R/R\y)^ the optimal values of R and Ry/ diverge and 
local charge neutrality is recovered in the mixed phase, where the energy density equals that of 
the Maxwell construction (see Fig. 19.171) . 

9.4.5 Thermal and neutrino trapping effects 

In following we discuss the thermal and neutrino trapping effects on the quark-hadron phase 
transition. In the left panel of Fig. 19. 18] we show the dependencies of the isothermal free-energy 
on the cell size for several temperatures in Eq. ()9.58p at ps = 2 po> using the constant bag 
constant and the surface tension as B = 100 MeV-fm~^ and T5 = 40 MeV-fm~^, respectively. 

Here we assume the non-uniform structure is the droplet and do not take into account neu- 
trinos to see the thermal effects. The quark volume fraction [R/Ry/)^ is fixed to be the optimal 
value at T = MeV for each curve. From the this panel, we can see that the minimum point is 
shifted to a larger value of R as temperature is increased, and eventually disappears for T > 60 
MeV. The temperature dependence of the free energy comes from the Coulomb energy, the 
surface energy and the correlation energy. By comparing the temperature dependence of these 
contributions, we can see that the correlation energy is primarily responsible to the behavior 
of the minimum point. It means that the mixed phase becomes less stable as temperature is 
increased [23] . 

After the minimization of the free energy, we can get the thermodynamic quantities, such 
as pressure, entropy, etc..@In the right panel of Fig. 19.181 we show the EOS (baryon density 
vs. pressure for the neutrino-free case at T = 30 MeV. Clearly, the pressure becomes to close 
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the droplet case at ps = 2 po ^ind different temperatures. The free energy is normalized by the 
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We also show, for comparison, the mixed phase by the Maxwell construction by thin solid line. 



to the one given by the Maxwell construction at finite temperature, since the mixed phase 
becomes unstable and its density regime is thereby largely restricted. In Fig. 19.191 we show the 
dependencies of neutrino trappings on the free energy of the mixed phase. Here, we adopt the 
density-dependent bag constant, and fixed the quark volume fraction {R/Riy)"^ at Y^^ = 0.01 
for each curve. For Yi,^ > 0.1, the minimum points disappear. The right panel shows the each 
component of the free energy. From this figure, we can see that both the correlation energy and 
the Coulomb energy mainly contribute to the behavior of the minimum point in the presence of 
neutrinos, while it is a only the correlation energy for the temperature dependence. 

To elucidate the neutrino-trapping effect in the mixed phase, we show the density profiles 
and the Coulomb potential for the slab case in Fig. 19.201 For comparison, we use the same 
cell-size here. High neutrino fraction enhances the electrons to satisfy the chemical equilibrium 
(see He — Pz/ in Eq. (I9.58P ). As a result, the Coulomb potential changes drastically at the 
trapped-neutrino case as shown in this figure. 

We do not present the pressure or EOS at the neutrino-trapped case as in the case of 
finite temperature because, at the realistic situations in proto-neutron stars as T =40 MeV, 
Yi{= Ye + Y^^) = 0.4, the mixed phase becomes unstable. Related to this discussion, the recent 
study suggested a new possible mixed phase for the neutrino trapped matter [36]. They imposed 
the local charge neutrality a priori instead of directly solving the Poisson equation, infered from 
our previous result in Fig. 19.191 In such case, the surface tension and the Coulomb energy do 
not work to make any pasta structure, and the mixed phase may resemble the amorphous phase 
of quarks and hadrons. Our results look to show an establishment of the local charge neutrality, 
but details will be discussed elsewhere. 
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they are shown per baryon. 



9.5 Summary and concluding remarks 

We have discussed general features of the phase transitions in compact stars, where multi- 
component substance are in equilibrium. These phase transitions are non-congruent, and the 
usual Maxwell construction cannot be applied any more. Using the bulk calculation, we have 
presented basic concepts of the non-congruent transition. Although these concepts are helpful 
in discussing and analyzing the properties of the phase transitions in compact stars, we must 
bear in mind that they are obtained in an idealized situation and must be carefully applied 
in realistic calculations by taking into account rather complicated effects such as the finite-size 
effects. Actually we found that uniform nuclear matter can exist in the mechanically unstable 
region. 

When the finite-size effects are taken into account, the structured mixed phase (pasta phase) 
emerges to form the crystalline structures. Charge screening effect and the rearrangement effect 
of particle densities in the presence of the Coulomb interaction is also emphasized; they some- 
times cause an instability of the geometrical structures in the mixed phase. In the extreme case 
the EOS resembles the one given by the Maxwell construction. 

Liquid-gas transition in asymmetric nuclear matter and hadron-quark deconfinement tran- 
sition are studied and their mixed phases have been calculated by using the Thomas-Fermi ap- 
proximation for particle densities and the Wigner-Seitz (WS) cell approximation for the pasta 
structure. We have shown many interesting aspects of these phase transitions. However, we 
have also seen that the maximum mass of hybrid stars is rather low in the light of a recent 
observation [73], since contamination of hyperons considerably softens EOS. The transition to 
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quark matter suppresses their appearance, but the maximum mass is still low due to the rather 
soft EOS of quark matter. So, we need another idea about EOS's of hadrons and quarks to 
circumvent this situation. 

So far, most works on the pasta structures have used the WS cell with ansatz about the 
geometrical structures like droplet, rod, slab and so on. If we want to see the possibility of 
the intermediate shapes, or study the mixed phase in a more realistic way, a study is desirable 
without the WS cell approximation or ansatz about the geometrical structures. Recently we 
have performed a three-dimensional calculation of non-uniform nuclear matter based on the 
relativistic mean-field model and the Thomas- Fermi approximation |50] . Introducing a cubic 
cell with periodic boundary conditions and discretizing it into grid, we numerically solve the 
coupled field equations for meson mean-fields and the Coulomb potential. Randomly distributing 
fermions (n, p, e) on the grid points as an initial condition, we relax their density distributions 
to attain the uniformity of the chemical potentials. As a preliminary result we demonstrate 
that pasta structures, previously given by the WS cell approximation, also appear in the three- 
dimensional calculation. 

Formation of nuclear pasta has been discussed in supernovae |7j or possible observation of 
the signals of the mixed phase has been suggested in the spectra of the gravitational waves 
[8]. The phenomenological implications of the pasta structures should deserve more elaborate 
studies; their thermodynamical or dynamical property is interesting; e.g. the neutrino opacity 
or the thermal conductivity inay be affected by the pasta structures, and their elasticity may 
cause a discontinuous change of internal structure of compact stars. In particular the crystalline 
structure in the core region may bring about novel dynamical effects on compact-star phenomena 
by way of deformation or oscillation. 
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